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DESIGN  OF  A  fflGH  POWER  X-BAND 
MAGNICON  AMPLIFIER 


I.  Introduction 

The  magnicon  [1'3],  a  “scanning  beam”  microwave  amplifier  related  to  the 
gyrocon  [4],  is  a  possible  replacement  for  klystron  amplifiers  in  future  electron  linear 
accelerators  operating  at  frequencies  in  the  X-band  and  above.  Compared  to  the 
klystron,  the  magnicon  offers  the  potential  of  higher  power  operation  at  frequencies 
greater  than  10  GHz,  higher  efficiency,  and  reduced  sensitivity  to  variations  in  load 
impedance,  while  preserving  the  traditional  klystron  attributes  of  high  gain  and 
frequency  and  phase  stability.  These  characteristics  may  also  make  the  magnicon  an 
attractive  device  at  lower  frequencies  (0.5-10  GHz),  for  a  variety  of  applications  in 
which  frequency  and  phase  stability,  combined  with  high  overall  efficiency,  are 
important  criteria.  Such  applications  include  high  average  power  proton  accelerators 
for  tritium  production  and  radioactive  waste  transmutation  (see  Ref.  [5]  and 
references  therein),  and  microwave  power  beaming  [6]. 

The  magnicon  consists  of  a  number  of  TM-mode  deflection  cavities  (a  drive 
cavity  and  a  series  of  gain  cavities),  followed  by  an  output  cavity.  The  deflection 
cavities  spin  up  an  electron  beam  to  high  transverse  momentum,  i.e.,  a  ■ 
where  v±  and  are  the  velocity  components  perpendicular  and  parallel  to  the  axial 
magnetic  field.  They  are  typically  designed  so  that  the  electron  transit  time  is  half  an 
if  period.  However,  the  final  deflection  cavity,  also  known  as  the  penultimate  cavity, 
may  consist  of  several  sections,  coupled  by  irises,  in  order  to  reduce  the  rf  fields 
required  to  spin  the  beam  up  to  high  a.  The  deflection  cavities  employ  a  rotating 
TMiio  mode,  producing  a  gyrating  electron  beam  whose  centroid  rotates  about  the 
Mtnuscripi  q>iiroved  November  16, 1993. 
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cavity  axis  at  the  drive  frequency.  In  addition,  the  sense  of  all  electron  velocities  and 
displacements  change  in  synchronism  with  the  advance  in  phase  of  the  rf  wave.  The 
output  cavity  extracts  principally  the  transverse  electron  momentum  in  a  gyroresonant 
“fast-wave”  interaction  that  extends  over  several  wavelengths.  It  employs  a  mode 
that  rotates  synchronously  with  the  beam  centroid  motion,  making  possible  a  highly 
efficient  output  interaction  without  requiring  beam  bunching. 

In  order  to  extrapolate  magnicon  operation  to  high  frequencies  (>5  GHz),  it  is 
possible  to  operate  the  output  cavity  at  a  multiple  of  the  frequency  of  the  deflection 
cavities  [2,  7,  81.  In  this  case,  the  output  cavity  will  employ  a  TM^io  mode,  where  m 
is  an  integer  greater  than  1.  This  permits  larger  cavities,  and  lower  rf  fields,  in  the 
deflection  system  than  would  be  required  if  the  entire  device  operated  at  the  output 
frequency.  However,  the  penalty  is  a  likely  lowering  of  the  efficiency,  due  to  the 
weaker  rf  fields  near  the  axis  of  the  output  cavity  for  m>l.  The  most  practical  value  of 
m  for  a  frequency  multiplying  magnicon  is  m-2.  In  this  case,  the  output  cavity 
operates  in  the  TM210  mode  at  twice  the  drive  frequency.  An  added  advantage  of  the 
frequency-doubling  configuration  is  that  approximately  the  same  magnetic  field  can  be 
employed  in  the  deflection  cavities  and  in  the  output  cavity.  This  is  because  the 
optimum  magnetic  field  in  the  deflection  cavities  conesponds  to  a  gyrofrequency 
approximately  twice  the  operating  frequency  of  the  deflection  cavities  [1],  while  in  the 
output  cavity  it  corresponds  to  a  gyrofrequency  iqiproximately  equal  to  the  operating 
frequency  [2].  A  generic  frequency  doubling  magnicon  configuration  is  illustrated  in 
Fig.  1. 

This  paper  presents  the  results  of  a  design  study  to  produce  a  circuit  for  an 
11.4  GHz  frequency-doubling  magnicon  amplifier  experiment  on  the  NRL  Long-Pulse 
Accelerator  (LPA)  Facility  19].  The  study  assumes  that  the  experiment  will  make 
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use  of  certain  existing  laboratory  equipment.  This  includes  a  single-shot  Marx 
generator  power  supply  driving  an  explosive  emission  diode  to  produce  a  -500  keV, 
172  A,  5.5-mm-diameter  solid  electron  beam  [10].  The  effect  of  finite  beam  diameter 
was  a  particular  concern  in  this  study.  In  addition,  the  study  assumes  that  the 
magnetic  field  will  be  produced  by  a  set  of  existing  pancake  magnets.  These  magnets 
are  3-cm  thick,  with  a  coil  inner  and  outer  radius  of  7.9  cm  and  15.75  cm.  They  are 
typically  separated  by  3. 3-cm  spacers.  However,  the  spacing  may  be  varied  to  tailor 
the  magnetic  field  profile.  Alternatively,  the  magnet  may  be  separated  into  several 
sections,  which  may  be  powered  by  separate  power  supplies. 

n.  THEORETICAL  MODELING  Al-ID  DESIGN  PROCEDURE 
A.  Physical  Model 

This  study  makes  use  of  steady-state  particle  simulation  employing  the 
realistic  fields  of  magnicon  cavities  connected  by  beam  tunnels,  and  includes  the 
effects  of  finite  electron  beam  diameter.  However,  it  does  not  model  either  space 
charge  or  beam  emittance.  Following  the  analytical  determination  of  approximate 
cavity  dimensions  and  operating  parameters,  the  simulation  procedure  begins  by 
solving  for  the  eigenmodes  of  particular  cavity  configurations,  including  the  effects  of 
beam  tunnels.  A  particle-pusher  code  then  integrates  the  trajectories  of  electrons 
through  these  fields,  first  of  individual  cavities,  and  then  of  an  entire  magnicon  device. 
The  particle  code  uses  a  steady-state  model  that  makes  use  of  the  phase  synchronism 
between  the  beam  and  the  wave.  Thus,  the  entire  simulation  is  carried  out  by 
following  a  single  temporal  slice  of  electrons,  one  macroparticle  deep  in  z. 

The  modeling  of  the  magnicon  must  take  into  account  the  modifications  to  the 
ideal  cavity  rf  modes  caused  by  the  presence  of  their  connecting  beam  tunnels.  These 
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“non-ideal”  ei^ects  can  strongly  influence  magnicon  operation,  and  their  evaluation 
requires  the  use  of  accurate  rf-fleld  solutions.  Figure  2  shows  rf  electric  ficJ  lines  in 
the  gain,  penultimate,  and  output  cavities  of  the  magnicon  circuit.  One  effect  is  a 
change  to  the  effeaive  interaction  length  of  each  cavity  and  to  the  rf-field-fcee  regions 
separating  the  cavities.  However,  this  can  be  compensated  for  in  the  design  of  the 
magnicon  cavities.  A  more  significant  effect  is  the  introduction  of  strong  transverse 
fields  and  radial  gradients  in  the  vicinity  of  the  cavity  beam  tunnel  and  irises.  These 
fields  and  gradients,  combined  with  a  finite  electron  beam  radius,  will  produce  velocity 
spreads  that  reduce  the  coherence  of  the  electron  motion.  These  effects  are 
particularly  important  in  the  penultimate  and  output  cavities,  where  the  combination  of 
initial  beam  diameter  and  high  induced  transverse  momentum  requires  the  largest 
beam  mnnels. 

Accurate  rf  fields  are  also  required  to  model  the  w-mode  of  the  penultimate 
cavity,  and  to  model  field  leakage  between  the  penultimate  and  output  cavities.  In 
order  to  account  for  this  last  effect,  the  penultimate  and  output  cavities  are  considered 
as  a  single  rf  structure,  which  is  solved  separately  for  the  dipole  operating  mode  of  the 
penultimate  cavity  and  the  quadrupole  operating  mode  of  the  output  cavity  (see  Fig. 
2b, c).  These  separate  solutions  are  then  superposed  in  the  particle  simulation. 

B.  Design  Goal 

The  objective  of  this  study  was  to  design  a  magnicon  circuit  with  maximum 
efficiency  for  the  parameters  of  the  proposed  NRL  experiment.  In  addition,  the  circuit 
should  provide  for  high  gain  (requiring  no  more  than  -1  kW  of  drive  power)  and  for  an 
effective  Q  no  higher  than  ~1(XX),  to  provide  a  reasonable  frequency  bandwidth  for  the 
device.  FinaUy,  the  circuit  should  be  stable  against  oscillation. 
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C.  Description  of  computer  models 

The  magnicon  codes  are  designed  to  calculate  realistic  if  fields  in  the  cavities 
and  to  determine  the  beam  electron  trajectories  in  those  fields. 

Calculation  of  the  fields  for  dipole  or  quadrupole  modes  of  axisymmetric 
cavities  with  arbitrary  shapes  is  performed  in  two  stages.  First,  two  components  of 
the  rf  field  are  calculated  using  the  MULTIMODE  code  [11].  In  this  code,  the 
eigenvalue  problem  is  solved  using  the  r  and  z  components  of  the  magnetic  field.  Hr 
and  In  this  case,  there  are  no  nonphysical  solutions.  Then  the  last  component, 
is  determined  from  the  equation 

+  (1) 

where  0)  is  the  eigenfiequency  and  c  is  the  speed  of  light.  The  term  (V^H)^  contains 
the  known  component  Hr .  This  procedure  does  not  decrease  the  precision  of  H^.  In 
the  MULTIMODE  code,  eight-node  second-order  isoparametric  fmite  elements  are 
used. 

Next,  it  is  necessary  to  determine  the  components  of  the  electric  field.  It  is 
possible  to  do  this  numerically  by  evaluating  the  derivatives  of  the  components  of  the 
magnetic  Held  or  by  solving  Maxwell’s  equations  using  a  fmite  element  method. 
However,  both  of  these  methods  give  unacceptably  large  errors.  Instead,  we  use  a 
boundary  integral  method  [12].  Using  the  MULTIMODE  code,  we  determine  the 
cavity  surface  current  density  distribution: 

j=nxH  (2) 
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where  n  is  the  unit  vector  normal  to  the  surface.  Using  the  current  density 
distribution,  we  can  find  the  vector  potential  at  any  point  r  inside  the  cavity: 


(3) 


where  k  =  m/c  ,  po  is  the  permeability  of  free  space,  7'  is  the  surface  source 
coordinate,  dS'  denotes  the  element  on  the  surface,  and  the  integral  is  over  the  entire 
cavity  surface.  Then,  it  is  possible  to  find  the  magnetic  field: 


W=/t^*VxA.  (4) 

Since  only  the  kernel  of  Eq.  (3)  depends  on  r,  it  is  possible  to  find  H  by  analytical 
differentiation  of  this  kernel.  The  electric  field  is  also  determined  from  Maxwell’s 
equations  by  analytical  differentiation: 

E=— VxVxA  (5) 

i0)£o 

Here,  £o  is  the  permittivity  of  free  space.  In  this  case,  all  six  field  components  will 
satisfy  Maxwell’s  equations,  and  have  the  same  precision  of  calculation.  For  the 
approximation  of  the  surface  current  density,  we  use  a  third-order  spline  fit.  To 
increase  the  speed  of  the  field  calculation  near  the  cavity  axis,  we  use  a  second-order 
paraxial  expansion  for  all  rf  field  components. 
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The  equations  of  motion  for  the  beam  electrons  are  integrated  using  the  third- 
order  Runga-Kutta  method.  During  the  integration,  energy  losses  and  cavity 
detunings  are  also  calculated.  A  finite  diameter  electron  beam  is  broken  up  into  a 
number  of  macroparticles,  spread  uniformly  over  the  beam  cross  section.  Thirty-seven 
macropaiticles  are  typically  employed,  but  up  to  177  macropaiticles  have  been  used  to 
verify  the  accuracy  of  the  final  simulations.  The  beam  may  be  modeled  either  with 
zero  initial  pitch  angle  (to  model  a  magnetic-field-immersed  diode)  or  zero  canonical 
angular  momentum  (to  model  a  shielded  diode). 

D.  Design  Procedure 

The  solution  for  the  realistic  cavity  rf  fields  is  an  input  to  the  particle  simulation 
code,  which  assumes  a  particular  amplitude  and  phase  of  the  fields  in  each  cavity,  and 
integrates  the  full  Lorentz  equations  for  the  electrons  through  these  fields.  Iteration  is 
used  to  find  self-consistent  values  of  amplitude  and  phase.  The  correct  rf  phase  in  a 
particular  cavity  is  assumed  to  be  the  phase  that  extracts  the  maximum  energy  from 
the  electrons.  For  the  gain  and  penultimate  cavities,  the  rf  amplitude  is  determined  by 
power  balance,  including  the  effect  of  a  finite  Q  due  to  waU  losses.  For  the  output 
cavity,  a  different  procedure  is  followed.  First,  an  rf  field  amplitude  is  assumed.  Next, 
the  rf  phase  is  adjusted  to  optimize  the  efficiency  for  the  chosen  amplitude.  Finally, 
the  output  Q  that  will  yield  overall  power  balance  is  calculated. 

The  process  begins  with  solutions  for  only  the  penultimate  and  output  cavities. 
These  preliminary  solutions  look  for  an  efficient  interaction  in  the  output  cavity  as  a 
function  of  the  amplitude  and  phase  of  the  rf  fields  of  both  the  output  and  penultimate 
cavities.  The  beam  entering  the  penultimate  cavity  is  assumed  to  be  unperturbed  by 
the  previous  cavities.  Separately,  the  gain  of  the  deflection  cavities  is  optimized  as  a 
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function  of  their  length  and  spacing.  Interesting  cases  are  then  made  self-consistent 
by  proceeding  backward  from  the  penultimate  cavity  toward  the  drive  cavity,  one 
cavity  at  a  time,  to  determine  the  field  amplitude  in  each  preceding  cavity  that  will 
drive  the  beam  in  order  to  produce  power  balance  at  the  fields  required  in  the  next 
cavity.  Finally,  to  ensure  self-consistency  and  overall  power  balance  through  the 
entire  device,  the  cavities  are  solved  again,  sequentially,  beginning  with  a  particular 
signal  in  the  drive  cavity. 

A  separate  area  of  concern  is  stability  against  self-excitation  (oscillation).  The 
principal  area  of  concern  is  the  penultimate  cavity.  The  magnicon  simulation  code  can 
be  used  to  study  the  stability  properties  of  magnicon  cavities  for  rotating  modes.  This 
is  done  by  looking  at  energy  balance  in  a  cavity  for  an  initially  unperturbed  electron 
beam  and  a  “small”  cavity  rf  field.  If  the  beam  loses  more  energy  to  the  cavity  mode 
than  the  mode  dissipates  in  wall  losses,  ibt  mode  is  unstable  and  will  grow  from 
noise.  The  same  energy  balance  calculation,  using  higher  cavity  rf  fields,  can 
determine  whether  the  oscillation  will  reach  large  amplitude,  or  whether  there  are 
unstable  states  of  “bard  excitation,”  that  might  be  excited  during  the  rise  of  the 
voltage  pulse.  However,  the  magnicon  simulation  code  can  only  model  synchronously 
rotating  modes,  and  thus  cannot  examine  the  TMqio  (klystron)  mode,  which  occurs  at 
approximately  3.6  GHz  in  the  deflection  cavities.  Symmetric  (m=0)  modes,  such  as 
the  TMqio  mode,  must  be  evaluated  by  another  means. 

In  order  to  estimate  the  penultimate  cavity’s  klystron-mode  instability 
threshold,  it  is  modeled  as  two  coupled  cavities.  The  problem  of  excitation  of  cavities 
coupled  by  small  holes  has  been  treated  by  Bethe  [13].  If  coi  and  eo2  represent 
eigenfrequencies  of  the  individual  cavities  of  the  penultimate  cavity,  the 
eigenfi’equencies  of  the  zero  and  n  modes  of  the  combined  system  are  given  by  o)  and 
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oKl+ic),  where  jc  is  the  coupling  coefficient  and  o)  If  ic«l6)i-<U2l/tu ,  the 

cavities  are  considered  to  be  weakly  coupled.  In  this  limit,  the  frequencies  of  the  two 
coupled-cavity  modes  will  be  approximately  equal  to  the  frequencies  of  the  two 
uncoupled  cavities.  In  this  case,  each  mode  will  have  one  filled  cavity  and  one  almost 
empty  cavity,  and  the  cavities  will  oscillate  independently.  In  this  small  jc  limit,  the 
threshold  for  oscillation  may  be  estimated  by  assuming  that  instability  can  occur  only 
if  the  cavity  detuning  caused  by  the  beam  is  of  the  same  order  as  the  difference  in  the 
cavity  eigenfrequencies. 

TMoiO“roo<le  fluctuations  in  the  first  cavity  will  modulate  the  beam  momentum, 
resulting  in  current  density  modulation  in  both  of  the  cavities.  This  modulation  causes 
both  frequency  detuning  and  a  change  in  the  beam  energy.  Specifically,  there  will  be 
electron  acceleration  in  the  first  cavity  (i.e.,  cavity  loading)  and  deceleration  in  the 
second  cavity  (cavity  excitation).  The  relative  detuning  A©/©  caused  by  the  beam 
may  be  estimated  (see  the  Appendix)  as: 

A©/©  -  R/Q  ■  Im(y), 

where  R  is  the  cavity  shunt  impedance,  Q  is  the  quality  factor,  and  Im(10  is  the 
susceptance  of  the  beam.  The  detuning  for  the  first  cavity  should  be  smaller  than  for 
the  second,  since  the  beam  is  not  yet  modulated.  To  estimate  the  instability 
threshold,  we  will  make  the  "worst  case”  assumption,  that  the  detuning  of  both 
cavities  will  be  of  the  same  order  and  of  opposite  sign.  The  magnitude  of  Im(10  cannot 
exceed  the  modulus  of  the  total  beam  admittance,  IT).  It  is  easy  to  show  (see  the 
Appendix)  that  for  a  relativistic  beam,  ll1»/o(fi)L/c)/()3^y^t/o).  Here,  L  is  the  length 
of  each  cavity,  Iq  is  the  beam  current,  where  mo  is  the  rest  mass  and  e  is 
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the  charge  of  an  electron,  y  is  the  relativistic  factor,  fie  is  the  electron  velocity,  and 
ffl/27t»3.6  GHz  is  the  rf  ftequency  of  the  TMqi  mode.  Then,  the  threshold  current  may 
be  estimated  as; 


Q  c 


(6) 


The  study  reported  in  this  paper  makes  use  only  of  steady-state,  or  time- 
independent  modeling.  Some  time-dependent  issues  of  magnicon  modeling  have  been 
explored  elsewhere  [14]. 


in.  DISCUSSION  OF  THE  Final  Circuit 

The  final  magnicon  circuit  includes  a  drive  cavity,  two  half-wavelength  gain 
cavities,  and  a  two-section  n-mode  penultimate  cavity,  all  operating  at  5.7  GHz,  and 
fmally  an  output  cavity  operating  at  11.4  GHz.  This  contiguration  is  illustrated 
schematically  in  Fig.  3.  Figure  4  shows  the  same  cavities  to  scale,  drawn  as  an 
overlapping  sequence  of  the  discrete  rf  structures  depicted  in  Fig.  2.  Figure  4  also 
shows  the  position  of  the  12  pancake  magnets.  These  are  arranged  in  two  sections, 
each  containing  6  magnets  separated  by  3.3-cm  spacers.  The  two  sections  are 
separated  by  a  larger  6.6-cm  spacer.  The  upper  part  of  Fig.  4  shows  the  resulting 
magnetic  field  profile.  The  design  parameters  are  summarized  in  Table  I. 

Figure  5  illustrates  the  “ideal”  operation  of  the  magnicon  by  analyzing  the 
case  of  a  single,  initially  on-axis  electron.  Notice  that  the  electron  transverse 
momentum,  plotted  as  the  “pitch  angle”  or  tan-^a,  progressively  increases  through 
the  deflection  cavities,  but  that  most  of  the  increase  occurs  in  the  penultimate  cavity. 
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Finally,  the  transverse  momentum  is  reduced  to  near  zero  in  the  output  cavity  and  at 
the  same  time  the  electron  energy  decreases  by  -58%.  (For  a  single  electron, 
substantially  higher  efficiencies  are  possible  if  a  is  allowed  to  be  »1.) 

Particular  areas  of  the  circuit  are  discussed  in  separate  sections  below. 

A.  The  Drive  and  Gain  Cavities 

The  spin  up  of  the  electron  beam  takes  place  in  a  series  of  deflection  cavities 
that  progressively  amplify  the  rf  signal  injected  into  the  drive  cavity  (see  Figs.  5  and 
6).  The  drive  and  gain  cavities  are  simple  cylindrical  cavities  utilizing  a  rotating 
TMno  mode.  The  calculated  ohmic  Q  of  these  cavities  is  ~15(XX).  The  drive  cavity 
will  be  loaded  by  external  coupling  pins,  in  order  to  match  the  drive  signal  into  the 
beam-loaded  cavity,  while  the  gain  cavities  are  loaded  principally  by  the  beam.  The 
use  of  beam  loading  to  lower  the  effective  Q  is  discussed  in  Section  IVB.  The  optimum 
cavity  length  corresponds  to  an  electron  transit  time  of  half  an  rf  period  (n  radians). 
The  diameter  of  the  drift  tubes  is  1  cm  in  this  region  of  the  magnicon  circuit.  The 
magnetic  field  was  chosen  to  be  1.6  times  the  cyclotron  resonant  value  at  500  keV. 
The  simulation  calculates  that  the  beam  will  shift  the  cold-cavity  resonant  frequency 
by  -15  MHz,  and  lower  the  Q  to  -1000.  The  calculated  transverse  electric  and 
magnetic  fields  along  the  axis  of  the  second  gain  cavity  are  shown  in  Fig.  7.  The 
transverse  electric  fields  are  the  result  of  the  beam  mnnels  at  either  end  of  the  cavity. 

B.  The  Penultimate  Cavity 

The  increase  of  rf  fields  is  approximately  exponential  in  the  cavities  of  the 
deflection  system  (see  Fig.  6).  As  a  result,  the  highest  fields,  and  most  of  the  beam 
spin-up,  will  occur  in  the  last  of  these  cavities,  the  penultimate  cavity.  If  this  cavity  is 
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similar  in  design  to  the  other  deflection  cavities,  very  high  rf  fields  are  required  to 
reach  a  beam  a  -I.  This  increases  the  likelihood  of  microwave  breakdown,  and  also 
can  lead  to  significant  wall  losses,  which  will  lower  the  overall  electronic  efficiency  of 
the  magnicon. 

In  order  to  lower  the  required  fields,  it  is  desirable  to  have  an  interaction  that 
can  be  sustained  in  a  favorable  rf  phase  for  more  than  half  of  an  if  period.  In  particular, 
the  use  of  a  multisection  cavity,  in  which  the  fields  change  by  n  after  each  half-rf- 
period,  will  serve  this  purpose  [2,  7].  In  principle,  a  two-section  cavity  should  need 
only  half  of  the  peak  rf  fields,  reducing  the  likelihood  of  breakdown.  Furthermore,  since 
wall  heating  is  proportional  to  the  square  of  the  rf  fields,  such  a  two-section  cavity 
would  dissipate  approximately  half  as  much  rf  power  as  a  single-section  cavity. 

A  further  advants^e  in  reducing  the  rf  fields  and  the  wall  losses  would  be 
gained  by  using  more  than  two  sections.  However,  a  stability  constraint  militates 
against  this.  Studies  have  shown  that  three-section  cavities  become  unstable  to 
oscillation  in  the  TMqiq  ic/2  mode,  in  which  the  three  cavities  act  approximately  as  a 
buncher,  drift  space,  and  catcher,  in  a  manner  analogous  to  a  klystron  [8].  Additional 
sections  would  only  make  matters  worse. 

Even  a  two-section  cavity  is  significantly  more  likely  to  be  unstable  than  a 
single-section  cavity.  The  stability  of  the  penultimate  cavity  against  oscillation  in  the 
TMi  10  n-mode  requires  a  lower  magnetic  field  than  the  optimum  field  in  either  the  half¬ 
wavelength  gain  cavities  or  the  output  cavity.  In  order  to  stabilize  the  penultimate 
cavity,  the  magnet  was  divided  into  two  sections  separated  by  6.6  cm.  The 
penultimate  cavity  is  located  in  the  vicinity  of  the  magnetic  field  minimum  produced  by 
this  separation  (see  Fig.  4). 
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The  stability  of  the  penultimate  cavity  against  oscillation  in  the  TMqio  nuxle  is 
estimated  using  Eq.  (6).  The  SUPERLANS  code  [IS]  is  used  to  calculate  the 
properties  of  the  TMoi  modes  of  the  penultimate  cavity.  For  our  case,  the  calculated 
value  of  RIQ  is  70  Q.  Also.  mi/2ic  =  3626  MHz  and  Q^/2n  =  3705  MHz.  In  addition, 
Vs  0.251  mi-n)2l/m,  implying  weak  coupling  as  discussed  in  Section  IID.  Thus,  by 
Eq.  (6),  the  threshold  current  is  /q  ~  440  A,  more  than  twice  the  design  value  of  the 
current.  One  should  note  that  a  similar  penultimate  cavity  has  been  used  for  the 
7  GHz  magnicon  [7,8]  at  the  Budker  Institute  of  Nuclear  Physics  (INP),  where  the 
beam  current  is  240  A  and  the  estimated  threshold  current  is  400  A.  In  this  device, 
the  klystron-type  instability  has  not  been  observed  [8]. 

The  two  sections  of  the  penultimate  cavity  are  of  equal  length  in  the  final 
circuit.  (Cavities  with  unequal  lengths  were  also  evaluated.)  The  cavity  is  designed 
so  that  each  of  the  sections  will  yield  approximately  the  same  absolute  value  for  the 
integral  of  the  transverse  magnetic  Held  along  the  cavity  axis,  i.e.,  a  **balanced** 
cavity.  This  is  accomplished  by  employing  a  slight  difference  (-0.5%)  in  the  cavity 
inner  radius  of  the  two  sections,  to  correct  for  the  asymmetry  in  the  beam  tunnel 
diameters  at  either  end  of  the  cavity.  The  transverse  electric  and  magnetic  fields  of 
the  penultimate  cavity  mode  as  a  function  of  position  along  the  axis  is  shown  in  Fig.  8. 
The  simulations  predict  a  maximum  if  electric  field  of  250  kV/cm  on  the  walls  of  the 
penultimate  cavity.  They  also  predict  that  the  beam  will  detune  the  cold  cavity 
resonance  by  -10  MHz. 

C.  The  Output  Cavity 

The  output  cavity  of  a  frequency-doubling  magnicon  employs  a  quadrupole 
TM210  mode,  which  rotates  synchronously  with  the  TMhq  mode  of  the  deflection 
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system.  The  cavity  length  is  chosen  as  a  compromise  between  the  high  fields  needed 
for  an  efficient  interaction  in  short  cavities,  and  the  lower  efficiency  that  results  from 
velocity  spreads  (and  other  effects)  in  longer  cavities.  The  design  calls  for  a  length  of 
5  cm.  This  corresponds  to  approximately  4  cyclotron  periods  of  the  electrons, 
assuming  an  initial  a-1  and  ~  0.8,  where  £2c  is  the  relativistic  cyclotron 
frequency  and  m  is  the  angular  frequency  of  the  radiation.  The  maximum  rf  electric 
fields  are  -250  kV/cm. 

Because  of  the  scanning  beam  synchronism  upon  which  the  magnicon  is  based, 
the  electron- wave  interaction  is  invariant  through  an  rf  period.  Thus,  for  an  ideal 
beam,  perfea  phase  synchronism  is  obtained  without  the  requirement  for  any  bunching 
mechanism.  However,  despite  the  invariance  in  time,  there  is  no  guarantee  that  all 
the  electrons  emitted  at  a  particular  time  will  remain  well  bunched  in  position  and 
phase  as  they  traverse  the  magnicon  ciicuit,  in  order  to  produce  an  efficient  interaction 
in  the  output  cavity.  The  latter  requirement  is  ultimately  a  constraint  on  the  initial 
beam  quality  and  on  the  spatial  gradients  of  the  cavity  rf  fields  over  the  cross  section 
of  the  elearon  beam.  In  addition,  space  charge  may  degrade  the  efficiency. 

Figures  9  and  10  show  optimized  simulations  of  the  final  magnicon  circuit  for  a 
S.S-rom-diameter  “magnetized”  and  a  2-mm-diameter  “shielded”  beam.  In  Fig.  9,  the 
interaction  of  a  5.5-mm-diameter  beam  with  the  cavity  fields  results  in  a  large  spread 
in  the  electron  trajectories.  A  significant  energy  spread  is  first  seen  in  the  second  gain 
cavity.  However,  the  most  significant  phase  mixing  occurs  in  the  penultimate  cavity, 
where  the  electrons  develop  a  large  spread  in  energy  and  in  pitch  angle.  The  nonideal 
efiects  in  both  of  these  cavities  result  in  large  part  from  the  large  radial  electric  fields 
shown  in  Figs.  7  and  8.  As  a  result,  the  interaction  in  the  output  cavity  is  relatively 
incoherent,  and  the  final  efficiency  is  only  24%,  much  less  than  58%  seen  in  a 
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comparable  single  particle  simulation  (see  Fig.  5).  In  Fig.  10,  for  the  2-mm-diameter 
shielded  beam,  the  entire  interaction  is  much  more  coherent,  resembling  the  single- 
particle  simulation  of  Fig.  S.  As  a  result,  the  final  efficiency  is  56%. 

The  2-mm-diameter  beam  corresponds  to  the  approximate  parameters  of  the 
INP  magnicon  gun  [2,  16].  The  correct  modeling  of  the  interaction  of  a  near-BnIlouin 
beam  such  as  that  produced  by  the  INP  gun  cannot  be  carried  out  without  including 
space  charge,  a  capability  that  does  not  presently  exist.  The  best  model  to  use  in  the 
absence  of  space  charge  is  not  clear.  However,  a  comparable  simulation  for  a  2-mm- 
diameter  magnetized  beam  yields  only  a  slightly  lower  52%  efficiency. 

Complete  simulations  were  also  carried  out  for  a  1. 2-mm-diameter  shielded 
beam  and  a  4-mm-diameter  magnetized  beam.  Figure  11  plots  the  efficiency  as  a 
function  of  beam  diameter  for  all  these  simulations.  It  is  noteworthy  that  the  efficiency 
falls  by  only  a  few  percent  up  to  beam  diameters  of  2  mm,  but  decreases  rapidly  for 
larger  diameters.  Figure  12  shows  the  efficiency  of  the  output  cavity  as  a  function  of 
the  rf  fields  in  the  penultimate  cavity  for  the  case  of  a  5  J-mm-diameter  beam.  The 
magnicon  efficiency  is  apparently  insensitive  to  small  variations  of  the  penultimate 
cavity  fields. 

IV.  technical  Problems  AND  DESIGN  Features 
A.  Rf  Fields  and  Breakdown 

The  axial  electric  field  of  the  TMno  mode,  in  combination  with  the  strong  axial 
magnetic  field  required  for  the  electron-wave  interaction,  creates  conditions  favorable 
for  multipactor  breakdown  in  the  deflection  cavities.  However,  this  problem  is  no 
different  in  principle  from  that  of  many  other  devices,  such  as  klystrons.  Nevertheless, 
it  may  be  an  area  of  concern  for  single-shot,  short-pulse  experiments  employing 
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relatively  poor  (10"5  Torr)  vacuum  pressures.  Effects  attributed  to  multipactor  were 
observed  in  a  preliminary  experiment  that  measured  the  gain  between  two  deflection 
cavities  [10].  A  second  area  of  concern  is  true  high-Held  microwave  breakdown  (Oeld 
emission).  In  order  to  avoid  breakdown,  the  peak  field  levels  in  both  the  penultimate 
cavity  and  the  output  cavity  are  well  below  the  Kilpatrick  limit  [17]. 

B.  Beam  Loading  (Gyrotropic  Effect) 

The  interaction  between  the  (initially  axicentered)  electron  beam  and  the 
TM 1 10  mode  of  the  deflection  cavities  is  highly  sensitive  to  the  value  of  the  axial 
magnetic  Held,  which  determines  the  sign  and  amplitude  of  the  energy  exchange 
between  the  electrons  and  the  mode.  In  addition,  the  sense  of  the  axial  magnetic  field 
breaks  the  symmetry  of  the  interaction,  resulting  in  different  behavior  for  the  two 
orthogonal  circular  polarizations  [1,2,18].  The  “preferred”  circular  polarization,  which 
produces  high  efficiency,  co-rotates  with  the  sense  of  the  electron  gyration. 

The  beam  loading  of  the  drive  cavity,  for  the  preferred  circular  polarization,  is 
identically  zero  for  an  on-axis  electron  with  no  initial  transverse  momentum  when  the 
relativistic  electron  cyclotron  frequency  in  the  axial  magnetic  field  is  twice  the 
operating  rf  frequency.  In  this  case,  the  electron  trajectory  remains  in  the  £^=0  plane 
of  the  rotating  rf  mode.  The  loading  is  also  near  zero  for  a  finite  electron  beam.  As  a 
result,  the  electrons  can  be  deflected  in  the  fint  cavity  without  gaining  or  losing 
energy.  However,  the  beam  loading  of  the  drive  and  gain  cavities  progressively 
increases  as  the  magnetic  field  is  lowered  from  this  value.  This  provides  a  valuable 
means  of  controlling  the  effective  cavity  Q  without  requiring  the  use  of  external 
coupling  or  lossy  materials.  Another  effect  of  beam  loading  is  to  alter  the  resonant 
frequency  of  the  cavities,  as  the  reactive  part  of  the  beam-circuit  interaction  down- 
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(up-)  shifts  the  resonance  from  the  vacuum  value  for  the  preferred  (opposite) 
polarization  [1,10].  This  shift  must  be  accounted  for  in  the  cavity  design. 

C.  Energy  Extraction 

The  Q-iiciot  of  the  output  cavity  is  determined  principally  by  the  output 
coupling.  The  output  coupling  is  chosen  to  yield  the  desired  electric  fields,  in  power 
balance  with  the  beam,  to  achieve  the  optimum  efficiency.  In  the  case  of  the  S.5-mm- 
diameter  beam,  this  requires  a  g  of  ~800,  while  for  the  2-mm-diameter  beam,  Q  should 
be  ~3(X).  In  comparison,  the  calculated  ohmic  0  of  a  copper  output  cavity  is  ~20,0(X). 
The  output  coupling  may  be  achieved  in  several  different  geometries.  The  principal 
variants  are  side  coupling  and  end  coupling. 

Side  coupling  is  the  method  conventionally  used  in  klystrons.  Its  advantages 
are:  1)  the  if  is  separated  from  the  beam,  simplifying  collector  design,  and  2)  the 
power  may  be  extracted  directly  to  in-band  rectangular  waveguide,  making  it 
immediately  available  for  applications.  There  must  be  at  least  two  side  couplers, 
separated  by  either  45®  or  135®  degrees,  to  provide  equal  loading  to  the  two 
orthogonal  linearly  polarized  components  of  the  circularly  polarized  TM210  mode.  The 
principal  disadvantage  of  side  coupling  is  the  potential  to  break  the  symmetry  of  the 
quadmpole  mode,  causing  conversion  to  other  modes.  Mode  conversion  into  dipole 
modes  will  lead  to  undesirable  leakage  of  microwave  power  into  the  penultimate 
cavity.  To  minimize  this,  additional  wall  perturbations  are  needed,  at  90®  intervals 
from  each  of  the  couplers,  to  preserve  quadmpole  synunetry. 

The  alternative  is  to  use  end  coupling.  For  end  coupling,  the  downstream  end- 
wall  of  the  output  cavity  will  contain  an  iris,  followed  by  a  circular  output  waveguide, 
rather  than  a  cut-off  beam  tunnel.  This  preserves  the  quadmpole  symmetry.  The  iris 
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diameter  and  thickness  are  chosen  to  achieve  the  desired  cavity  value.  However, 
the  microwave  mode  must  pass  through  the  collector  before  reaching  the  output 
window,  and  will  also  need  further  processing  to  produce  a  useful  mode,  e.g.,  for 
coupling  to  an  accelerating  structure.  This  is  the  same  situation  faced  by  most 
gyrotron  devices,  including  gyroklystrons  [19],  although  for  those  devices  TE  rather 
than  TM  modes  are  involved.  The  use  of  end  coupling  is  envisioned  for  the  first  test  of 
the  magnicon  circuit. 

V.  Summary 

Using  a  magnicon  simulation  code,  we  have  designed  a  microwave  circuit  for  a 
high  power  frequency-doubling  magnicon  amplifier  at  1 1.4  GHz.  The  design  study  was 
carried  out  using  the  approximate  parameters  (500  keV,  172  A,  5.5-mm-diameter)  of 
an  electron  beam  produced  from  a  field-emission  diode  on  the  NRL  LPA  Facility. 
Using  these  parameters,  the  predicted  magnicon  efficiency  is  24%.  Based  on  a  large 
number  of  simulations,  this  design  appears  to  be  insensitive  to  minor  variations  of  the 
various  experimental  parameters.  However,  the  final  efficiency  depends  strongly  on 
the  properties  of  the  electron  beam,  with  a  large  premium  to  be  gained  from  reducing 
the  initial  beam  diameter.  The  same  circuit,  if  driven  by  a  beam  produced  by  the  INP 
magnicon  gun,  should  produce  an  efficiency  exceeding  50%.  The  initial  experiment  will 
use  the  existing  single-shot  LPA  beam.  However,  further  development  of  the 
magnicon  as  a  practical  source  in  this  frequency  range  will  require  the  use  of  a  rep¬ 
rated  modulator  combined  with  a  high  quality  electron  beam  produced  by  a  thermionic 
diode. 
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Appendix 


The  purpose  of  this  appendix  is  to  outline  the  derivation  of  the  expression  for 
the  relativistic  beam  admittance  employed  in  Section  II .  The  admittance  represents 
the  beam  loading  on  the  if  circuit.  This  derivation  follows  along  the  lines  indicated  in 
[20].  To  determine  the  several  contributions  to  the  admittance,  we  shall  first  obtain 
the  current  modulation  in  the  second  part  of  the  penultimate  cavity  (hereafter  referred 
to  as  the  second  cavity)  resulting  from  the  velocity  modulation  in  the  first  part  of  the 
penultimate  cavity  (hereafter  referred  to  as  the  first  cavity).  For  simplicity,  the  motion 
is  assumed  to  be  one-dimensional,  and  the  cavities  are  represented  by  short  gaps,  of 
length  d,  the  centers  of  which  ate  separated  by  a  distance  L. 

The  relativistic  equation  of  motion  is 

(A1> 


where  7  is  the  relativistic  mass  factor,  mg  is  the  electron  rest  mass,  7  is  the  electron 
velocity  along  the  z  axis,  and  V  is  the  potential  difference  across  the  first  cavity.  From 
the  definition  of  y-  (1  -  it  follows  that 


(A2) 


Combining  Eqs.  (Al)  and  (A2),  one  obtains 


(A3) 
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For  the  in  the  first  cavity, 


^V»h.san<ot  ,  (A4) 

az  d 

where  Vi  is  the  peak  voltage  across  the  gap.  Making  use  of  Eqs.  (A3)  and  (A4),  and 
following  along  the  lines  indicated  in  [20],  the  ac  current  induced  in  the  second  cavity 
is  given,  in  phasor  notation,  by: 


where 


is  the  beam-coupling  coefficient.  In  phasor  notation,  the  voltage  across  the  first  g^  is 
-jV\,  and  hence  the  electronic  admittance  is  given  by  Ji2 


where /5s=t?/c.  In  our  case,  mL/n-tc. 
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Following  a  similar  procedure,  the  small-signal  admittance,  l^i,  of  an 
unmodulated  beam  traversing  the  distance  d  across  the  gap  in  the  Hrst  cavity  is 
readily  determined.  The  result  may  be  expressed  in  the  form 


2[l-cos(tu<//t;)]-(m<f/i;  )sm(mrf/z7 ) 

Re(^i)  =  So - V/  si - 

2{(o  div) 


2sin(Q)<//p  )|l+cos(Q)rf/p  )] 


2(fi)df/p  ) 


where 


«0  = 


2/n 


Note  that  for  Af-l,  \Yi2\  and  IFil  are  of  the  same  order  of  magnitude. 
The  total  beam  admittance,  Y,  may  be  written  in  the  form 

r.o*(cc-±). 


in  terms  of  the  lumped  parameters  of  the  equivalent  parallel  RLC  circuit,  with  G=R-^. 
The  beam  loading  imparts  a  frequency  detuning  Ao>  to  the  circuit  with  respect  to  the 
resonant  frequency,  to  =  (LC)“*^,  given  by 


Aq) 

o> 
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where  Q  is  the  cold  cavity  quality  factor.  Since  the  two  parts  of  the  penultimate  cavity 
are  assumed  to  impart  detunings  of  the  same  order  of  magnitude,  but  of  opposite 
signs,  the  net  detuning  is  given  by  the  expression  indicated  in  Section  II . 

To  assess  the  effect  of  space  charge  on  the  electron  bunching  due  to  a  beam  of 
finite  radius,  we  solve  the  one-dimensional  relativistic  Vlasov  equation  for  the 
electron  distribution  /: 

^  +  o^-kl^^  =  0  (A5) 

at  ai  op 

where  p=yniQV  is  the  axial  momentum  and  E  is  the  axial  electric  field.  Expanding 
about  the  equilibrium  distribution,  Eq.  (A5)  relates  the  perturbed  current,  /,  to  the 
electric  field.  In  the  limit  in  which  the  wavelength  is  large  compared  to  bfy,  where  b  is 
the  cavity  radius,  we  can  also  relate  the  electric  field  in  the  cavity  to  the  current  via: 

where  rb  is  the  electron  beam  radius  and  Zo^jiQ/eQ)^^  is  the  impedance  of  free  space. 
Assuming  that  the  perturbation  varies  as  txpijm-jkz),  the  effect  of  space-charge 
waves  on  the  beam  is  represented  by  the  wavenumber 
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where  U Q=mQc'^l\e\.  For  the  5.5-rmn-diameter  beam  and  cavity  radius  b=3.2  cm, 
ikpxO.Ob  ale,  indicating  that  space  charge  has  a  minor  effect  on  the  bunching. 
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Table  I 


Magnicon  Design  Parameters 

Input;  -1  kW  O  5.7  GHz 

Beam:  500  keV,  172  A 

Deflection  cavities;  3  half- wavelength  TM  no-mode  cavities  (1  drive,  2  gain) 

Cavity  RadiuS'-3.2  cm;  Cavity  Length=2.3  cm; 

Beam  Tunnel  Radius=0.5  cm;  Beam  Tunnel  Length=1.9  cm 

Penultimate  cavity:  n-mode  two-section  TMno  cavity, 

Cavity  Radius~3.2  cm.  Cavity  Length=5  cm 

Output  cavity:  TM210  mode,  Q-800  (for  5.5-mm-diam.  beam) 

Radius~2.15  cm;  Lengths5  cm 

Maximum  if  fields  at  walls:  -250  kV/cm 

(Penultimate  and 

Output  Cavities) 

Output:  fout  =  4/slaC  =  1 1  GHz 

P  -  20  MW  @  ■n~24%  for  5.5-mm-diam.  beam 
P  -  50  MW  @  11-56%  for  2-mm-diam.  beam 
Overall  gain  -  44-48  dB 
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Figure  1,  Schematic  of  frequency-doubling  magnicon  configuration. 
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Figure  4.  Schematic  of  the  cavities  imd  magnetic  field  coils,  including  the  axial 
magnetic  field  proHle  used  in  the  simulations. 
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Figure  7.  Transverse  electric  and  ms^netic  fields  on  the  axis  of  the  second  gain 
cavity  as  a  function  of  position  along  the  axis.  The  cavity  outline  is  also 
shown. 
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Figure  8.  Transverse  electric  and  magnetic  fields  on  the  axis  of  the  penultimate 
cavity  mode  as  a  function  of  position  along  the  axis. 


37 


Radius  (cm)  Pitch  Angle  (degrees)  Electron  Energy  (keV) 


800 


600 


400 


200 

0 

60 


40 


20 

0 

3.2 


1.6 


0.0 

40.00  46.52  53.04  59.56  66.08  72.60 

z  (cm) 

Figure  9.  Simulation  of  the  flnal  magnicon  design  for  a  5.5-mm-diameter  beam, 
showing  the  spatial  evolution  of  electron  energy  and  velocity  pitch  angle, 
and  the  electron  trajectory. 
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Figure  10.  Simulation  of  the  final  magnicon  design  for  a  2-mm-diameter  beam  from  a 
shielded  cathode. 
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Figure  11.  Magnicon  efficiency  versus  beam  diameter.  Shielded  cathode  (•)  and 
magnetized  beams  (  ■  )  cases  arc  included.  Note  that  the  point  at  0  mm 
is  the  same  for  shielded  and  magnetized  beams. 
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Figure  12.  Efficiency  of  the  output  cavity  versus  field  in  the  penultimate  cavity. 


